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Suffice it to remember what Kant asserted; 
that progress in every science is measured in 
terms of its use of mathematics. (Gori 2004, 44) 


1. Introduction 


International relations investigates a vast universe of political phenomena, 
most of it constituted by a mix of continuity and discreteness. The duration of 
diplomatic relations among countries, of peace between states, of international 
treaties, and of global international regimes in diverse policy domains are con- 
tinuous variables; as are distance between capitals, speed of great power tran- 
sitions, and probabilities associated with all international events. By contrast, 
the formal composition of a country’s diplomatic organization, of alliances, 
governmental and nongovernmental international organizations, as well as the 
requisites of effective deterrence and other extant policies, are discrete varia- 
bles. Time, space, territories, and emotions are generally continuous, but with 
discrete features such as barriers, thresholds, empty spaces, layers, and bounda- 
ries, which are discrete. This hybrid texture of continuity and discreteness—i.e., 
“concreteness,” meaning simultaneously continuous and discrete—is ubiquitous, 
consequential, and fundamental in international relations, as reflected by theory 
and research across the discipline (and throughout social science in general). 

Research and analysis of hybrid theories of international relations is conduct- 
ed through mathematical tools from the infinitesimal calculus of Newton and 
Leibniz and discrete calculus developed in recent decades. Both are needed to 
understand real-world international phenomena that are otherwise not know- 
able through purely historical or narrative discussions (Gillespie 1976; Kline 
1985; Cioffi 1998). Until recently, however, infinitesimal and discrete calculi 
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have remained largely disjoint. Here we demonstrate a new unified calculus of 
hybrid functions with novel applications to a small, albeit representative and 
convincing sample of international relations theories. As a scientific system for 
exploration and discovery, this new analysis uncovers novel, significant, and 
often surprising features and properties of international phenomena that are 
otherwise inaccessible and, therefore, remain unknown, through earlier ap- 
proaches. This investigation demonstrates results through formal mathemati- 
cal and computational analysis supported by visual analytics, similar to the use 
of alternative “diagnostic imagery” in medical analyses or different ensembles 
of observational instruments in scientific research. 

The next section provides examples of hybrid phenomena in international 
relations, followed by a section on the methodology of nabladot calculus for 
unified hybrid analysis. The fourth section investigates three specific cases that 
demonstrate hybrid analysis applied to international phenomena. Since our inter- 
est is substantive (asin all applied mathematics), we focus mainly on significant 
features of international political phenomena rather than purely mathematical 
themes. The last section provides concluding remarks. 


2. Hybridity and hybrid functions in international relations theories 


Continuity and discreteness—ontological hybridity—are present in the 
following international phenomena and their respective theoretical explanans: 

Peace and other compound international events. All international events— 
e.g, political integration, alliance formation, conditions for peace, success of 
international regimes, deterrence requirements, nuclear proliferation contain- 
ment, and international communication—are compound events, in the sense of 
probability theory, because they are always caused by several (i.e., more than 
one) conjunctive events (Bittinger and Crown 1982; Bruschi 1990; Goertz and 
Starr 2003). Consequently, the probability ofan international event is a function 
of some discrete number of causal conditions required for its occurrence and a 
continuous value of probability associated with each causal event (Wohlstetter 
1968; Cioffi 1998, chs. 5-7). 

Growth of great powers or empires (Taagepera’s law). At the actor-level of 
analysis, as an empire expands from some initially small size up to its maxi- 
mum size, its growth is governed by a logistic function, where time is strictly 
continuous but the rate of polity expansion is discrete—by chunks of territory 
(provinces, other administrative units, conquered territories bound by natural 
barriers, fortifications, and other limiting factors (Taagepera 1968; 1978; 1979). 

Wright-Snyder crisis theory of war. An inter-state war never “comes out of 
the blue” but, rather, originates from a prior crisis (a crisis being a metastable 
phase-transition, in complexity-theoretic terminology). The crisis- and bargain- 
ing-based theory explains the onset of interstate war as the violent escalation 
outcome ofa process with multiple outcomes, war amongst them. In this theory, 
the probability of war in a given epoch is a function of some discrete number of 
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crises during the epoch and the probability of escalation to war in each crisis 
(Wright 1942, 1271-76, fn. 38; Deutsch 1978; Snyder and Diesing 1977, 13-5 
et passim; Cioffi 1998, 160-63). 

Loss-of-power gradient. At the relational or meso-level of analysis, continu- 
ity and discreteness over accessible time and space have mixed effects on the ex- 
ercise of power (“power projection”) at distances away from home base. While 
the rate of decay can assume positive continuous values, distance from home 
base is discrete (determined by military bases, supply chains, and other discrete 
systems and networks), so the overall function of these two quantities is hybrid 
(Boulding 1962; Wohlstetter 1968). 

Size of war alliances (Horvath-Foster law). The frequency of war alliances 
in politico-military history shows a pattern that decreases with the size of the 
alliance. Informally, there have been many small war alliances, very few large 
ones, and an intermediate number in between (Horvath & Foster 1963). This 
is known as a discrete Yule-Simon distribution with continuous parameter and 
is symptomatic of complex systems and generative processes that are far from 
equilibrium; otherwise, the size distribution of war alliances would be normal 
or Gaussian (as in the height or weight of persons). 

Warfare and international systemic polarity (Midlarsky’s law). At the sys- 
temic or macro-level of analysis, the annual frequency of warfare in the interna- 
tional system varies in proportion to the number of great powers in the system, 
known as polarity. However, the frequency of wars increases with marginally 
decreasing increments in systemic polarity, so this too is a hybrid, nonlinear re- 
lationship (Midlarsky 1974). 

Numerous other instances of international phenomena and corresponding 
theoretical explanations exist in international relations. Here we shall use in- 
stances 1, 3, and 6 to demonstrate how and why they are investigated through 
nabladot analysis, as described in the next section. 


3. Mathematical methods for hybrid functions in IR 


Consider a hybrid function, 2 = y(X,Y),suchthat(X,Y) > Z € R,where 
Х and Ү are real-valued continuous and discrete independent variables, respec- 
tively.' Nabladot analysis of a hybrid function Ф(Х, У) begins by (1) clarify- 
ing the hybrid domain of y, specifically its substantive subdomains x € X 
and y € Y along each independent variable—such subdomain always being a 
bounded subspace of some broader mathematical domain and then, (2) specify- 
ing each variable’s unit of measurement. This initial phase of analysis normally 
includes various graphs of Ф for visual analysis, which are typically 2D or 3D 


We shall restrict attention to scalar hybrid functions, although vector hybrid functions also 
arise in nabladot analysis of scalar functions, as we shall see later in section 4, just as they do 
in classical analysis. Surveys of mathematical methods in political science and international 


relations include Cioffi (1979), Ashford et al. (1993), Moore and Siegel (2013). 
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surface graphs and contour plots of the hybrid function under investigation. As 
we shall see, ensembles of these interrelated graphs constitute theoretical land- 
scapes—complete with singularities, basins, escarpments, canyons, and other 
topographic features—that provide sometimes surprisingly faceted or nuanced 
explanations and deeper understanding of each of the hybrid functions. 

The next phase—and first properly analytical step in theoretical analysis—is 
to closely examine the causal effect of each independent variable on the depend- 
ent variable of interest, which is how the emergent field (dependent variable Z) 
is generated by the hybrid domain—given that y maps the former (causes) onto 
the latter (effects). This consists of two steps that examine absolute and stand- 
ardized effects, respectively. First, the first-order derivative and first-order dif- 
ference of hybrid function are separately calculated, graphed, and examined, to 
understand absolute variations with respect to changes in X (continuous inde- 
pendent variable) and Y (discrete). This phase maintains the original units of 
measurement corresponding to each variable, since derivatives and differences 
are simple rates of change. 

Second, the point elasticity and the arc elasticity of ©, denoted by nz and ny, 
respectively, are calculated to understand how patterns of variation in percentage 
change in each independent variable compare independent of units of measurement 
(which is what elasticity operators 7.(2)апа 7y(Z) are designed to investigate). 
Additional graphs and visual analytics are used as well to better understand the 
structure and effects of elasticities—and add to the theoretical landscape of each 
hybrid function. This second phase results in transformed standardized dimen- 
sional space without units of measurements, making all independent variables 
and their direct effect on the dependent variable directly comparable. These re- 
sults lead to one or more dominance principles, which are law-like statements 
that specify which independent variable has greatest causal effect on the depend- 
ent variable of the hybrid function—a fundamental property not always obvious 
from simple inspection of the hybrid function under investigation. 

The analytical process thus far has focused on scalar properties of the IR hy- 
brid function under investigation. The first nabladot operation is to calculate the 
hybrid gradient of Z to discover the magnitude and direction of changes in Z as 
a function of changes in X and Y. The result of applying the nabladot operator (a 
vector operator) to scalar hybrid function ¢ isa hybrid vector function Ф = Vy 
with x- and y-components. The hybrid gradient in two dimensions is the scalar 
vector product calculated using the new nabladot vector operator y (note the dot 
within the nabla symbol), which is defined as follows: 


Ур = 0,9 i+ Ayy j, (1) 


* Economists call this “comparative statics,” a phrase we shall not use here because time can 
be an independent variable of interest (e.g., as in Taagepera’s law of empires) which—by 
definition—is not static. 
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where, by convention, i andj denote unit vectors along x- and y-dimensions, re- 
spectively, and 0, and A, denote the first-order derivative and first-order diffe- 
rence with respect to X and Y (Cioffi 2014; 2017; 2019; 2020; 2020).° Note that 
the resulting nabladot gradient of hybrid function ø is a striated vector field 
with a first-order partial derivative component along the x-axis (continuous) 
and a first-order partial difference component along the y-axis (discrete), hence 
the striation of the vector field’s topology.* 

The absolute and standardized norms of the hybrid gradient Vy аге calcu- 
lated next, along with corresponding graphs for investigating the resulting vec- 
tor field. Each pair of plots for a vector field and corresponding norm should 
use identical domains to facilitate understanding through comparative analy- 
sis. Cardinal directions (N, E, $, W) are used for simple orientation in graphs. 
Other hybrid operations of nabladot calculus equivalent to the divergence, curl, 
Laplacian, Hessian, and Jacobian are subsequently calculated to shed addition- 
al (and usually new) light on the original function Z = y(X,Y).through the 
medium of nabladot operators, each supported by additional graphic analyses. 

The main results of nabladot analysis shed new light on fundamental, real- 
world, substantive properties and features of the original hybrid function under 
investigation, features that remain hidden or inaccessible through other forms 
of analysis. Each main formal expression is accompanied by an interpretation 
in plain English, although this is not always possible without some loss of preci- 
sion or clarity. Some results can be somewhat complicated nonlinear functions 
that do not further simplify; we prefer them that way rather than introducing 
artificial approximations which may be simpler but unrealistic or unnatural ob- 
jects, unlike real IR phenomena. In most cases an ensemble of images and visual 
analytics (Thomas and Cook 2005; Wellin 2013) of complicated functions can 
add significant clarity. 

Among the most important substantive (and testable) results from nabladot 
analysis are the dominance principles mentioned above—they explain which 
independent variable has dominant effect on the dependent variable, a major 
theoretical (and arguably policy) question impossible to answer ex ante—as 
well as other characteristic phenomena of interest (e.g., discrete striations, in- 
flection or “tipping” points, asymptotes and other singularities, constant or 
invariant subfields, and others) revealed by geometric and topological infor- 


3 ІК scholars rarely consider the presence of vectors in international relations, other than 


metaphorically. This analysis demonstrates the rigorous analysis of vectors and vector fields 
in IR using formal methods from nabladot calculus, as in the next section. To contain nota- 
tion, we shall use i and j to denote unit vectors along continuous and discrete dimensions, 
respectively, rather than create new unit vectors for each variable. 

Use of the partial derivative with respect to Y (a discrete variable) instead of the partial dif- 
ference—which is often used in approximations—produces a measurable error that varies 
in magnitude depending on the structure of 2 and values of Y. Measurable discrepancies 
between the two operators (nabladot and classical nabla) are demonstrable but beyond the 
present scope due to space limitations (Cioffi 2021). 
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mation. In addition, interesting scalar and vector fields of become accessible 
to direct investigation through formal tools of nabladot calculus and analysis. 
Anovel and valuable feature of this approach is that nabladot calculus provides 
exact results in analytical investigations where the classical infinitesimal calcu- 
lus of hybrid IR functions would provide approximations with errors over the 
discrete domain of independent variable(s). 


4. Applications to areas of international relations 


Here we shall investigate three illustrative cases of IR theories and research 
areas—numbered 1, 3, and 6in section 2—where hybrid functions play a central 
role in describing and explaining political phenomena. The scientific purpose is 
to deepen our understanding and provide foundations for more advanced analy- 
sis. Each “case study” follows the analytical procedure just outlined in section 3. 


4.1 Case 1: Peace and other international events 


All international events in the real world are “compound” because they are 
always produced bynecessary conditions specific to the event. Such causal neces- 
sity is universal—a fundamental axiom in all domains of international relations 
theory. For example, consider the event defined by the following expression: 


S = “a state of stable peace exists between two countries.” (2) 


This is a compound event because S requires the following set of causally 
necessary conditions, each of which constitutes an event by itself: 

1. Neither state will pursue issues deemed as highly threatening to the other, 
as opposed to acting completely oblivious or independently of extant fo- 
reign interests. 

2. Each state may prefer to negotiate over colliding interests, before escala- 
ting to war. 

3. When states do negotiate, they may—depending on conditions—find a 
nonviolent resolution. 


The universal existence of such necessary conditions makes S a compound 
event, by definition. 

Specifically, every international event E is produced by causal conjunction 
(operator Л) of necessary events—i.e., set-theoretic intersection (operator П) or 
Boolean logic product (operator AND). Let {%:}/_, denote a set of N necessary 
events that produce E. Causal production of E is specified by an event function, 
We : {Xi} > E, which maps necessary events in {%;} onto E using causal con- 
junctions. The number of events in a compound event (its “size” or “conjunctiv- 
ity”), is called cardinality, |E| = {1, 2,3,- -- ,n} C №, which is always a natural 
number (positive integer), so cardinality (or “event size”) is always a discrete 
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variable. For example, compound event 5 in equation 2 (“a size 3 event”) is first- 
order conjunctive with respect to its three causally necessary events.* Formally, 
we can summarize these ideas through the following expression: 


Ee P(E) = E A E2 A A En = A E (3) 


Next, an international event E has probability Pr( E) that is determined by 
the naturally uncertain occurrence of necessary conditions—a type of causa- 
tion known as probabilistic causality (Salmon 1980; Suppes 1984; Eels 1991). 
Based on Kolmogorov’s (1933) fundamental theorem of compound events, the 
joint probability E of compound event E is given by the following expression: 


Pr(E) = Pr(E,)- Pr(E2)-...- Pr( n) = Про, (4) 


where events E; are independent; when they are not independent, conditional 
probabilities are used and, by Kolmogorov’s theorem, they still multiply. Letting 
Pr(E) = E and Pr(E;) = pi, since all probabilities are continuous variables (i.e., 
Eand p are variables, not events), we can rewrite equation (4) in simpler notation: 


E=py-po-...-Pn =] [pi (5) 
i=1 
=P", (6) 
where Р is the probability of causal events and N = |E] is the cardinality or 


event size. Note that equation 6 is a hybrid function, because probability is con- 
tinuous оп the unit interval of real numbers [0,1]—which lies between (and in- 
cludes) impossibility (E = 0) and certainty (E = 1.0)—while event size is always 
a discrete or natural number (of necessary conditions). 

The politically relevant domain ofinterest—in this case—isboundby0 < P< 1 
and 2 < N 5 20 along causal probability and size dimensions in the hybrid (р, 
n) domain. Specifically, equation 6 is a hybrid function linearly dependent in P 
and exponentially in N, so this means that changes in either variables will cause 
different political effects on the probability of an event. 

The 3D surface graph of equation 6 is in Figure 1а, which shows values of 
the probability E of an international event rising through a steep escarpment 
as causal probability P increases and cardinality N decreases. The contour plot 


$ An event function Y is also called indicator function, structure function, or production func- 


tion, in areas of mathematics, engineering, and economics, respectively. 
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in Figure 1b shows the graph looking straight down, which highlights the am- 
ple basin floor where causal probability is very low (0< Р 0.9), flanked by the 
steep north-south escarpment along the east edge as Р—>1.0. 

Note that change caused by each variable differs, as shown in Figures 1с 
through f. Calculating the partial derivative and partial difference of E with re- 
spect to P and N, respectively: 


ðE = NPN-1 for P’s effect on Е, in Figures 1c and d) (7 
р 
А.Е = Ptt- PN (for N’s effect on E, in Figures le and f). (8) 


We see immediately that changes in causal probability P and event size N 
have different effects on the probability of an international event E, besides the 
trivial observation that any change on either variable has some effect on E (since 

дЬЕ + 0 and Ay E = 0, per equations 7-8). To wit: 

Opposite political effects. Whereas change in causal probability P has pro- 
portional or positive effect on E, change in event size N has an opposite effect. 
This is shown by the purple deep bottom in Figure 1с versus the high red plate- 
au in Figure le; both features observed over approximately the same underlying 
(p, n)-domain. 

First- and second-order effects. Whereas the contour plot of E(P, N) shows 
strictly concave isocontours (Figure 1b), the contours of the derivative 0,E 
(Figure 1d) also showa convex, mild spur protruding on the southwestern wall 
of the escarpment, near p = 0.75 (third contour, between green and blue eleva- 
tions), which means more complex change for small-size events. This is a se- 
cond-order effect and not at all intuitive from the basic model (equation 6 and 
Figures la and b). 

Second-order effects caused by changes in event size. The graphs of the chan- 
ge in E with respect to event size, An E (equation 8 and Figures le and f), show 
a pronounced ravine or canyon along the north-south direction dropping into a 
deep precipice at relatively high values of Pas N — 2. Interestingly, in this case 
all the isolines have mixed concavity (low P) and convexity (higher P values), 
which is another second-order effect. 

Geometrically opposite extrema. The extreme high range of д„Е in Figures 
1с and d and the extreme low range of АЕ in Figures le and f have opposite 
(or inverse) geometries with some common similarities: the former rises from 
a flat basin to an escarpment while the latter drops from a large plateau toward 
a deepening canyon that dives into a deep well. Both features are indicative of 
major political effects on event probability E caused by underlying changes in 
causal probabilities P and event size N. 
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Figure 1. Probability of an international event E(P, N) as a hybrid function of the 
probability P of N necessary conditions for its occurrence. (a) 3D surface of hybrid 
function E = PY; (b) contour plot of (а); (c) 3D surface of 0, Е, the first-order partial 
derivative of E(P, №) with respect to P; (d) contour plot of (c); (е) 3D surface of A, Ё, the 
first-order partial difference of E(P, N) with respect to N; (f) contour plot of (e). 


39 


‘These first insights begin to shed light on the nature of international events 
as a function of their necessary conditions. However, causal probability P and 
event size N are measured in different units (probability and number of events, 
respectively), so absolute rates of change (equations 7 and 8) fail to explain which 
variable has dominant effect on E. To solve this problem we obtain and analyze 
standardized (or percentage) rates of change in E with respect to Pand N (called 
elasticities), as in Figure 2. Calculating the point elasticity and the arc elasticity 
of E with respect to P and N, respectively: 


YE =N for P’s percentage effect on E, in Figures 1c and d 9 
р 
mE =(P—1)N (for N’s percentage effect on E, in Figures le and f). (10) 


We see immediately that these standardized effects on E are very different 
from the earlier absolute, unit-based effects (equations 7 and 8). 

The point elasticity—percentage change in the probability of an interna- 
tional event with respect to percentage change in P—is constant, as in Figures 
2a and b, meaning that a percentage variation in probability of an international 
event E relative to a percentage variation in causal probability P is determined 
solely by event size N and is independent of P. By contrast, arc elasticity—per- 
centage change in E with respect to percentage change in N—is linear in both 
independent variables, as in Figures 2c and d, resulting in a nonlinear scalar field 
as evidenced by the contour plot. 

Comparing the two elasticities (equations 9 and 10) answers the universal ques- 
tion concerning which of the two causal variables has the dominant or greater ef- 
fect on the probability ofan international event. Since N>(P-1)N, this means that 
point elasticity np is greater than the arc elasticity 7),,. Therefore, E is more sensi- 
tive to change in P (probability of necessary causal events) than to change in car- 
dinality N—a result that may be called dominance principle for international events. 

Although their effects differ, in reality both variables have joint, concurrent 
effects on international events. The joint effects of P and N on E can be better 
understood by calculating the gradient of E with respect to both variables using 
the nabladot operator, as follows: 


VE = 6,P%i+ An Pj (11) 
= NPN) i+ (P-1)PX j, (12) 


which is a two-dimensional vector function E = ¥(P, N). The emergent vector 
field of this hybrid gradient is in Figure 2e and corresponding vector magnitude 
or norm |E|(P, N), which is a scalar function, is shown in Figure 2f. We see im- 
mediately that norm |Е | is very similar to дЕ in Figures le and d, whichis not 
an intuitive result or obvious insight that could be obtained from casual com- 
parison between the two dissimilar equations for the arc elasticity (equation 8) 
and the magnitude of E. 

This concludes the first analysis of our three “case studies.” The next two are 
presented in slightly abbreviated form to omit some procedural repetitions while 
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maintaining the method outlined in section 3. Analysis ofthe probability ofinterna- 
tional events in this first case leads us to investigate the probability of war onset—a 
special class ofinternational event offundamental historical significance since early 
antiquity and of theoretical interest (atleast) since Thucydides—in the next section. 
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Figure 2. Elasticities and gradient of the probability ofan international event E. (a) 3D 
surface of point elasticity Np(E) respect to causal probability P; (b) contour plot of (a); 
(c) 3D surface of arc elasticity n(E) with respect to number of necessary conditions 
N; (d) contour plot of (с); (e) vector field of the dot-gradient vector function VE; (£) 
contour plot of (e). 
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4.2 Case 2: Crisis dynamics and onset of war 


Explaining the outbreak of international war as caused by а prior crisis esca- 
lation process—which includes challenging and resistance moves, bargaining, 
signaling, and other events, as opposed to some other causal mechanism—was 
first proposed by Quincy Wright (Wright 1942, 1271, fn. 38) and later extend- 
ed and generalized by Glenn Н. Snyder (Snyder Diesing 1977, 13-7 et passim). 
The frame of reference here is the inter-state relational level of analysis and the 
specific explanandum is the probability W ofa state being at war over a period of 
time, called an epoch in probability theory. During an epoch a country expe- 
riences a number C of crises (defined as episodes during which hostilities may 
occur), each individual crisis having probability w of escalating to war. As a re- 
sult, the probability of no war over C crises is (1 – w)° (by Kolmogorov’s theo- 
rem), so a country’s epochal probability of war is given by 


= 1 (1-а) (13) 


which is а bivariate nonlinear hybrid function, where w is continuous over ће 
closed unit probability interval [0, 1] and С> 2 is discrete. The case when C= 1 (a 
single crisis during an entire epoch) is trivial, since W(w, 1) = w, as is easily shown. 

The politically relevant domain is bound by 0 < w < 1 and 2 < C 5 20 (same 
as before, both equations being functions of compound events). 

The 3D surface graph of W(w, C) isin Figure 3a, which shows the probability 
of war rising rapidly to a maximal plateau as escalation probability w and C in- 
crease. The contour plot in Figure 3b looks straight down, which highlights the 
broad plateau where war probability converges to 1, flanked by the steep north- 
south escarpment along the west edge as w > 0. 

Each variable increases Win a different way, as shown in Figures 3c through 
f. Calculating the partial derivative and partial difference of W with respect to 
w and C, respectively: 


д =C(1—w)°"! (for w’s effect on W, in Figures Зс and d) (14) 
A.W =w(1—w)? (for C’s effect on W, in Figures Зе and f). (15) 


Changes in escalation probability w and number of crises C have clearly dif- 
ferent effects on epochal war probability W, although both functions are posi- 
tively valued: 

Congruent political effects. Change in either escalation probability w or num- 
ber of crises C has a direct effect on W, as shown by strictly positive values of the 
graphs of derivatives and differences in Figure 3c through f. 


6 


While Wright formulated the high-level crisis-probability framework, Snyder completed it 
by providing the probabilistic causal mechanism (Suppes 1984) within the crisis branching 
process. 
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First- and second-order effects. Whereas the contour plot of W(w, С) shows strictly 
concave isocontours (Figure 3b), contours of the derivative 3u W (Figure 3d) also 
show a convex, mild spur protruding on the southwestern wall of the escarpment, 
near p =0.2, a second-order effect in opposite direction to the previous case (Figure 
1d) and, again, not apparent from the basic model (equation 13 and Figures За апа). 
Second-order effects caused by changes in number of crises. The graphs of the 
change in E with respect to number of crises, A.W (equation 14 and Figures Зе and 
f), show a pronounced spur along the north-south direction descending from a high 
value of W along low values of w as C increases away from the minimal value of 2. All 
the isolines showmixed concavity and convexity, which is another second-order effect. 
Geometrically opposite spurs. The mild C-grown spur shown by 0,,W in 
Figures 3c and d and the more pronounced w-grown spur in A.W in Figures 
3e and f have a type of reflective symmetry: the former protrudes from the C 
boundary at lowvalues of w and C, while the latter extends from the w-boundary 
from low values (p ~ 0.2) toward high values of C. Both features are indicative 
of nonlinear effects on war probability W caused by underlying changes in cri- 
sis escalation w and number of crises C. 

Having obtained some initial insights on the nature of war probability as a 
function of crisis dynamics, we now investigate the Wright-Snyder hybrid model 
through standardized variables not based on units and, therefore, enable direct 
comparisons of causal effects. In this case we shall proceed by obtaining and 
analyzing elasticities of W with respect to w and C, as shown in Figure 4. Cal- 
culation of the point elasticity and the arc elasticity of W with respect to w and 
C yields the following set of hybrid equations: 


= C-1 › W in Fi 

NW = Сш(1- ш) “эша effect оп › ш Figures (16) 
1—-(1-w)? aan 

= 1 (for C’s percentage effect on W, in Figures (17 

Me Ca ane ) 4c and d) ( ) 


Based on these equations we see again that these standardized effects on W 
are quite different from the absolute, unit-based effects uncovered earlier (equa- 
tions 14 and 15 and associated figures). Recalling the meaning of elasticities, 
here, point elasticity stands for percentage change in the probability of war on- 
set with respect to percentage change in crisis escalation probability w, while 
arc elasticity is the percentage change in epochal war probability Wwith respect 
to percentage change in number of crises C during the epoch. In this case both 
elasticities are rather complicated rational hybrid functions, including denom- 
inators that are exponential in C (from the standardizing transformation), as 
shown in Figures 4a through d. While formal analysis is feasible, visual analyt- 
ics of graphs reveal numerous interesting features. 

Figures 4a and b both show that point elasticity is high (i.e., war probability 
is strongly affected by escalation probability) at low w values and highest at low- 
est values of both w and C (red levels of the escarpment). This means that the 
risk of war changes most when crises are few and escalation probability low. By 
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Figure 3. Epochal probability of war W(w, С) аз а hybrid function of crisis escalation 
probability w and number of crises C: computational imagery from visualization 
analytics. (a) 3D surface of the hybrid function W = (1 – w); (b) contour plot of (a); 
(c) 3D surface of O,,W, the first-order partial derivative of W(w, C) with respect to w; 
(d) contour plot of (с). (е) 3D surface of a „ the first-order partial difference of W(w, С) 
with respect to C; (f) contour plot of (e). 
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contrast, at the blue-green levels, war probability W is less sensitive to such in- 
stabilities (greater number of crises and higher escalation probability). 

A rather subtle and surprising feature occurs along the front edge of the 
3D surface in Figure 4, where point elasticity for all values of w at the minimal 
boundary of C = 2 is convex (bulging up), whereas elsewhere (away from the 
front edge of the surface) point elasticity is strictly concave (д) < 0). This makes 
minimal epochs with only two crises rather special, which is not immediately 
intuitive. This particular property vanishes in all epochs with a multiplicity of 
crises beyond just two (C > 3). Concavity in point elasticity of W with respect 
to w accelerates as the number of crises surpasses the first single digits. 

By contrast, arc elasticity of war probability Wis strictly concave, as seen in 
Figures 4c and d, convexity in this case being nonexistent (even at low C levels 
along the front edge). 

Comparing the two elasticities (equations 16 and 17) yields the following 
dominance principle: epochal probability of war W is more sensitive to change 
in escalation probability w than to change in number of crises C, because point 
elasticity y is greater than arc elasticity n . The different but joint effects of es- 
calation probability w and number of crises C on epochal war probability W are 
analyzed and understood by calculating the gradient of W with respect to both 
variables using the nabladot operator, as follows: 


VW = д, [1—(1—w)°]it+ Ac [1- (1-) 9] j (18) 
= С(1- ш)С-11+0(1- 0) j, (19) 


which is а two-dimensional vector function W = Ф (w, С). The resulting vector 
field of this hybrid gradient is in Figure 4e and corresponding vector magnitude 
or norm |W|(w, С), now a scalar function, is in Figure 4f. 

We see from these results that both probability vector fields and norms |W| 
and |E|, in Figures 2e and f and Figures 4e and f, resemble each other—a sur- 
prising result from comparative analysis—based on perfect bilateral vertical 
symmetries around the P = w = 0.5 axis. This is another political property not 
apparent from simple inspection of the basic models but consistent with formal 
fundamental symmetry and equivalence between causal logic conjunction and 
disjunction (De Morgan’s laws) associated with the probability of international 
events (AND-based conjunctive) and epochal war probability (OR-based dis- 
junctive), respectively. 


4.3 Case 3: Frequency of war and systemic polarity (Midlarsky’s law) 


Having just examined the probability of war at the relational level, in this 
last case study we again change our frame of reference, this time turning to the 
systemic level of analysis. The formal theoretical (and empirically supported) 
explanation for the annual frequency of warfare as a result of fundamental po- 
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Figure 4. Elasticities and gradient of the epochal probability of war W. (a) 3D surface 
of point elasticity n (W) with respect to crisis escalation probability w; (b) contour 
plot of (a); (c) 3D surface of arc elasticity y (W) with respect to epochal number of 
crises C; (d) contour plot of (c) ; (e) vector field of the dot-gradient vector function 


VW; 


46 


(f) contour plot of (e). 


litical uncertainty among great powers was pioneered by Мапи$ 1. Midlarsky in 
1974. This causal theory is based on aspects of strategic uncertainty rooted in 
Shannon’s information entropy and related concepts. The specific explanandum 
is the annual amount of warfare experienced in the international system, which 
is an emergent property generated not just by crises (asin the previous case) but 
also outright conquest, revenge, colonization, and all other types of wars. In any 
given year, Midlarsky’s theory predicted that annual global war frequency ¢ is 
determined by systemic polarity © among the great powers.’ Formally, 


¢=Klog®, (20) 


which is a hybrid bivariate nonlinear function, where K is a continuous propor- 
tionality parameter over the closed interval [1, 10], historically, and systemic po- 
larity © > 1 is discrete. The case when Ө = 1 (a single hegemonic power) is the 
most peaceful, since ф(Ө = 1) = 0, and the function's (К, @)-domain once again 
lies within Cartesian quadrant I. 

The 3D surface graph of equation 20 is illustrated in Figure 5a, which shows 
war frequency rising оп hilly slope as polarity © and K increase. The associated 
contour plot in Figure Sb looks “straight down hill,” highlighting the relatively 
mild gradient on the skirting slopes of the hybrid function, which is a surface 
bound by a linear slope along K but a logarithmically convex slope along Ө, a 
radial pattern observable in the plot’s isocontours. 

Each independent variable increases ф in a different way, as shown in Fig- 
ures Sa and b. Calculating the partial derivative and partial difference of ф with 
respect to Капа Ө, respectively: 


Ong = Ө (for K’s effect on W, in Figure 5c) (21) 
Лоф= K [In(O+1)—In®] (for O’s effect on ¢, in Figure Бе). (22) 


Here, too, we see that changes in parameter Капа polarity © have clearly dif- 
ferent effects on war frequency ф,1п this case through a simple univariate discrete 
function thatis strictly discrete in ©, and through а more complicated bivariate hy- 
brid function (Figures Sc and d) with two terms in both dimensions, respectively. 
Congruent political effects. Changes in parameter K or polarity © have pro- 
portional effects on ¢, as shown by strictly positive values of both graphs of de- 
rivatives and differences in Figures Sc (2D) and e (3D), respectively. The graph 
of 0;,¢ lacks a contour plot, since it has univariate domain in Ө. 
Isomorphism of ф and 0.0. War frequency ¢ and its rate of change with re- 
spect to Kare isomorphic, as shown by Figure Sc and the estimated hybrid fun- 


Midlarsky’s theory remains one of the most complete and empirically validated formal the- 
ories in international relations—and, surprisingly, one of the least known among conflict 
researchers. The initial validation based on 1815-1945 war onset data should not be retested 
on earlier and later datasets. 
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Figure 5. Annual frequency of war ф(К, ©) asa hybrid logarithmic function of systemic 
polarity Ө (Midlarsky’s law): (a) 3D surface of the hybrid function ф = K In Ө; 
(b) contour plot of (a); (с) 2D graph of 0;,¢, the first-order partial derivative of ф with 
respect to K; (d) plot of the original function fitted on historical data constrained 
by the logarithmic condition $(1) = 0; (e) 3D surface of A a the first-order partial 
difference of ф with respect to ©; (f) contour plot of (е). 
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ction ф(К,Ө) fitted to original data in Figure 54 (redrawn from Midlarsky 
1974, 420, fig. 3). 
Deceptive simplicity. Aside from the relatively simple convexity of war fre- 
quency ¢ induced by systemic polarity ©, this theoretical model seems rather 
uncomplicated, as would appear from diagnostic images in Figure 5. However, 
such an impression may be superficial, as shown next by the second part of the 
analysis on elasticities and vector fields. 

Finally, we can now investigate Midlarsky’s law using standardized variables, 
as in the previous cases. First, we calculate and analyze percentage change in ф 
with respect to Капа Ө, as shown by the elasticities in Figures ба and b. Calcu- 
lation of point elasticity and arc elasticity of ф with respect to parameter K and 
polarity © yields the following set of hybrid equations: 


кф = 1 (for К?з percentage effect on ф, in Figure ба) (23) 


йе “as 


mo 1 (for O’s percentage effect on ¢, in Figure 6b). (24) 

We immediately see that standardized effects on ф are again quite different 
from earlier unit-based results (cf. equations 21 and 22 and associated figures). 
Here, point elasticity represents percentage change in annual war frequency 
with respect to percentage change in parameter K, whereas arc elasticity meas- 
ures percentage change in ¢ with respect to percentage annual change in polar- 
ity ©. The former has a constant value of 1 while the latter is a rational hybrid 
function with logarithms of polarity, as shown in Figures 6a and b, respectively. 

Interestingly, arc elasticity (Figure 6b) exhibits a singularity under unipo- 
larity, where $(k, 1) = œ, indicating a major transition from unipolarity to bipo- 
larity. In this case arc elasticity continues a rapid drop with increasing polarity, 
whichis the systemic trend experienced in contemporary history of major pow- 
ers after the Soviet-American Cold War around 1989. 

Comparing the two elasticities (equations 23 and 24; cf. also their respective 
graphs in Figures 6a and b) yields the following dominance principle: annual 
frequency of war ¢ is more sensitive to percentage change in parameter K than 
to change in systemic polarity © under unipolar and bipolar systemic structures, 
but the reverse is true under tripolarity and higher-order structures—because 
the relationship reverses between © = 2 and 3. This is a surprising qualitative 
transition that is invisible in the original model but is clear once the dimensions 
are standardized by elasticities. 

The different and joint effects of parameter K and polarity © on war frequen- 
су ф can be seen by calculating the gradient of ф with respect to both variables 
using the nabladot operator, as follows: 


уф = & (KIn©)i+ Ag (KIn®)j (25) 
= log 0 i+ K [((nO0+1)—-Im9] j, (26) 
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which is a two-dimensional vector function Ф = ¥ (К, Ө). The resulting vector 
field of this hybrid gradient is in Figure 6c, which shows: (i) heterogeneity along 
the two dimensions; (ii) general southwest-northeast orientation; (iii) increas- 
ing intensity as K increases for low values of polarity with a hot spot in the NW 
corner; (iv) divergence from a line at approximately 60 degrees; and (5) curling 
associated with divergence. 

The corresponding vector magnitude or norm |Ф |,а scalar function, is shown 
by the contour plot in Figure 6d, on the same domain as the field. This shows 
other clear patterns, including a better view of (i) the hybrid gradient field Ф that 
drops from the NW region into the uni- and bi-polar basin at minimal values of 
K—clearly where the most peaceful worlds are found—and (ii) a distinct view 
of the high ridge beyond К> 5 and O > 2. 

Given such marked differences between elasticity functions, it is best to in- 
vestigate the standardized gradient, as follows: 


a ‹ ө | 
КЕТЕ (Ө + 1) А 
=1i+0 | ШӨ 1 j, (28) 


which is a two-dimensional vector function Ф* = ¥* (K, ©). The resulting vec- 
tor field of this hybrid standardized gradient is seen in Figure 6e, which shows 
differences that are best highlighted by the corresponding 3D plot in Figure 
6e. Here we see that the standardized gradient with respect to polarity under- 
goes a precipitous decline from unipolarity to bipolarity, after which it tapers 
off much more gradually and this feature is independent of K, consistent with 
elasticity results. 


5. Conclusions 


We began this chapter by observing that hybrid functions—formal models 
containing a combination of continuous and discrete variables—have been pre- 
sent in international relations theory since antiquity, “hiding in plain sight,” but 
their rigorous analysis has been impeded by the generally disjoint nature and 
established practices of infinitesimal calculus and discrete calculus. This prob- 
lem has now been resolved by a unified approach that is feasible and fruitful, as 
provided by nabladot calculus. 

Analysis of three separate cases in this chapter showed how and why, far 
from being intractable or only amenable to approximations, hybrid functions 
in international relations theory contain numerous interesting features and in- 
sightful properties that shed new light on our understanding of international 
phenomena. Nabladot analysis of each hybrid function—and subsequent com- 
parative analysis across them—in each case revealed previously unknown and 
often scientifically surprising theoretical landscapes ofinternational phenomena 
ranging from generic international events, to conditions of peace, crisis dynam- 
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Figure 6. Elasticities and gradients of the annual frequency of war ф. (a) point elasticity 
(Ф) with respect to scale parameter К; (b) arc elasticity (Ф) with respect to systemic 
polarity ©; (c) vector field of the dot-gradient vector function УФ; (d) contour plot of (с). 
(е) vector field of the standardized dot-gradient vector function Уф; (f) 3D plot of (e). 


ics, and warfare; all within the unified methods enabled by the hybrid nabladot 
operator and associated concepts, rather than through disjoint calculi or error- 
prone approximations. 

These new theoretical landscapes and research frontiers are exciting and 
their application is still in a preliminary but already a demonstrably promising 
stage. Dominance principles that rank the influence of causal independent vari- 
ables, singularities, previously undetected phase transitions, the deep nature of 
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probabilistic causality, and other scientifically intriguing results have intrinsic 
value for our understanding of international relations. As “progressive problem- 
shifts”, in the sense of Lakatos (1973; cf. also Gillespie 1976 and Moore 2001), 
this novel and emergent corpus of scientific knowledge also provides rich and 
creative foundations for more advanced analyses that enhance our theoretical 
as well as practical understanding. 
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